A convergent finite element method with adaptive √ 3 refinement

نویسندگان

  • Jan Maes
  • Adhemar Bultheel
چکیده

We develop an adaptive finite element method (AFEM) using piecewise linears on a sequence of triangulations obtained by adaptive √ 3 refinement. The motivation to consider √ 3 refinement stems from the fact that it is a slower topological refinement than the usual red or red-green refinement, and that it alternates the orientation of the refined triangles, such that certain features or singularities that are not aligned with the initial triangulation might be detected more quickly. On the other hand, the use of √ 3 refinement introduces the additional difficulty that the corresponding finite element spaces are nonnested. This makes the setting nonconforming. First we derive a BPX-type preconditioner for piecewise linears on the adaptively refined triangulations and we show that it gives rise to uniformly bounded condition numbers, so that we can solve the linear systems arising from the AFEM in an efficient way. Then we introduce the AFEM of Morin, Nochetto, and Siebert adapted to our special case for solving the Poisson equation. We prove that this adaptive strategy converges to the solution within any prescribed error tolerance in a finite number of steps. Finally we present some numerical experiments that show the optimality of both the BPX preconditioner and the AFEM. ∗Currently at Materialise, Leuven, Belgium †Department of Computer Science, K.U.Leuven, Celestijnenlaan 200A, 3001 Heverlee, Belgium A convergent finite element method with adaptive √ 3 refinement Jan Maes Adhemar Bultheel Department of Computer Science, Katholieke Universiteit Leuven, Belgium Abstract We develop an adaptive finite element method (AFEM) using piecewise linears on a sequence of triangulations obtained by adaptive √ 3 refinement. The motivation to consider √ 3 refinement stems from the fact that it is a slower topological refinement than the usual red or red-green refinement, and that it alternates the orientation of the refined triangles, such that certain features or singularities that are not aligned with the initial triangulation might be detected more quickly. On the other hand, the use of √ 3 refinement introduces the additional difficulty that the corresponding finite element spaces are nonnested. This makes the setting nonconforming. First we derive a BPX-type preconditioner for piecewise linears on the adaptively refined triangulations and we show that it gives rise to uniformly bounded condition numbers, so that we can solve the linear systems arising from the AFEM in an efficient way. Then we introduce the AFEM of Morin, Nochetto, and Siebert [21] adapted to our special case for solving the Poisson equation. We prove that this adaptive strategy converges to the solution within any prescribed error tolerance in a finite number of steps. Finally we present some numerical experiments that show the optimality of both the BPX preconditioner and the AFEM.We develop an adaptive finite element method (AFEM) using piecewise linears on a sequence of triangulations obtained by adaptive √ 3 refinement. The motivation to consider √ 3 refinement stems from the fact that it is a slower topological refinement than the usual red or red-green refinement, and that it alternates the orientation of the refined triangles, such that certain features or singularities that are not aligned with the initial triangulation might be detected more quickly. On the other hand, the use of √ 3 refinement introduces the additional difficulty that the corresponding finite element spaces are nonnested. This makes the setting nonconforming. First we derive a BPX-type preconditioner for piecewise linears on the adaptively refined triangulations and we show that it gives rise to uniformly bounded condition numbers, so that we can solve the linear systems arising from the AFEM in an efficient way. Then we introduce the AFEM of Morin, Nochetto, and Siebert [21] adapted to our special case for solving the Poisson equation. We prove that this adaptive strategy converges to the solution within any prescribed error tolerance in a finite number of steps. Finally we present some numerical experiments that show the optimality of both the BPX preconditioner and the AFEM.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Comparison of different numerical methods for calculating stress intensity factors in analysis of fractured structures

In this research, an efficient Galerkin Finite Volume Method (GFVM) along with the h–refinement adaptive process and post–processing error estimation analysis is presented for fracture analysis. The adaptive strategy is used to produce more accurate solution with the least computational cost. To investigate the accuracy and efficiency of the developed model, the GFVM is compared with two versio...

متن کامل

Stability and Convergence of Finite-Element Approximation Schemes for Harmonic Maps

Abstract. This article discusses stability and convergence of approximation schemes for harmonic maps. A finite element discretization of an iterative algorithm due to F. Alouges is introduced and shown to be stable and convergent in general only on acute type triangulations. An a posteriori criterion is proposed which allows to monitor sufficient conditions for weak convergence to a harmonic m...

متن کامل

A convergent adaptive stochastic Galerkin finite element method with quasi-optimal spatial meshes

We analyze a-posteriori error estimation and adaptive refinement algorithms for stochastic Galerkin Finite Element methods for countably-parametric, elliptic boundary value problems. A residual error estimator which separates the effects of gpc-Galerkin discretization in parameter space and of the Finite Element discretization in physical space in energy norm is established. It is proved that t...

متن کامل

Fast Finite Element Method Using Multi-Step Mesh Process

This paper introduces a new method for accelerating current sluggish FEM and improving memory demand in FEM problems with high node resolution or bulky structures. Like most of the numerical methods, FEM results to a matrix equation which normally has huge dimension. Breaking the main matrix equation into several smaller size matrices, the solving procedure can be accelerated. For implementing ...

متن کامل

A Nodal-based Finite Element Approximation of the Maxwell Problem Suitable for Singular Solutions

A new mixed finite element approximation of Maxwell’s problem is proposed, its main features being that it is based on a novel augmented formulation of the continuous problem and the introduction of a mesh dependent stabilizing term, which yields a very weak control on the divergence of the unknown. The method is shown to be stable and convergent in the natural H(curl; Ω) norm for this unknown....

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2009